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ABSTRACT 

If thick disks are ubiquitous and a natural product of disk galaxy formation and/or evolution 
processes, all undisturbed galaxies which have evolved during a significant fraction of a Hubble time 
should have a thick disk. The late-type spiral galaxy NGC 4244 has been reported as the only 
nearby edge-on galaxy without a confirmed thick disk. Using data from the Spitzer Survey of Stellar 
Structure in Galaxies (S'*G) we have identified signs of two disk components in this galaxy. The 
asymmetries between the light profiles on both sides of the mid-plane of NGC 4244 can be explained 
by a combination of the galaxy not being perfectly edge-on and a certain degree of opacity of the 
thin disk. We argue that the subtlety of the thick disk is a consequence of either a limited secular 
evolution in NGC 4244, a small fraction of stellar material in the fragments which built the galaxy, 
or a high amount of gaseous accretion after the formation of the galaxy. 

Subject headings: galaxies: individual (NGC 4244) — galaxies: photometry — galaxies: spiral — 
galaxies: structure 



1. INTRODUCTION 

Thick disks are detected in edge-on and very inclined 
galaxies as an excess of flux at some distance from the 
mid-plane, typically at a few thin disk scale-heights. 
They usually have exponential profiles with a larger 
scale-height than the 'canonical' or thin disks. Thick 
disks were first detected by Tsikoudi (1979) and were de- 
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fined by Burstein (1979). A thick disk has subsequently 
been found in the Milky Way (Gilmore & Reid 1983) and 
in basically all kinds of disk galaxies (see the introduction 
of Yoachim & Dalcanton 2008 and references therein) . 

Several theories of thick disk formation have been 
proposed. It has been suggested that thick disks 
are formed through mechanisms of internal evolution, 
namely disk heating by gravitational instabilities (Bour- 
naud et al. 2007; 2009) and/or stellar scattering by en- 
counters with features such as dark matter subhaloes, 
spiral arms or giant molecular clouds (Villumsen 1985; 
Hanninen & Flynn 2002; Benson ct al. 2004; Hayashi 
& Chiba 2006; Haywood 2008; Kazantzidis et al. 2008; 
Schonrich & Binney 2009a; Schonrich & Binney 2009b 
; Bournaud et al. 2009). In this picture the thick disk 
formation may take less than one gigayear in the case 
of fast dynamical heating by a very clumpy initial disk 
(Bournaud et al. 2009) or may be much slower (secu- 
lar evolution; Villumsen 1985; Hanninen & Flynn 2002; 
Benson et al. 2004; Kazantzidis et al. 2008). The specific 
type of secular evolution which causes a slow heating of 
the stellar orbits has been termed 'blurring' by Schonrich 
& Binney (2009a; 2009b). 

Another possibility is that the thick disk is formed 
through accretion of stellar material from merging satel- 
lites at the time of galaxy formation (Abadi et al. 2003; 
Yoachim & Dalcanton 2006). A third possibility is that 
the thick disk is a consequence of in situ star formation 
(Brook et al. 2004, Elmegreen & Elmegreen 2006) or of 
star formation with a high initial velocity dispersion in 
very massive star clusters (Kroupa 2002). 

Whatever the thick disk formation mechanism is, it 
has to act over most, if not all, disk galaxies as thick 
disks have been found to be ubiquitous (Dalcanton & 
Bernstein 2002; Yoachim & Dalcanton 2006). There- 
fore, if any disk galaxy was found not to possess a thick 
disk, it would be likely to have suffered a peculiar ori- 
gin and/or evolution. A search in the literature yields a 
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candidate, NGC 4244, which has been claimed either to 
be a doubtful case (Yoachim & Dalcanton 2006), or not 
to have a thick disk at ah (Fry et al. 1999). NGC 4244 
is a highly inclined apparently undisturbed Scd galaxy 
located 4.4 Mpc away (Seth et al. 2005a). At that dis- 
tance, one arcsecond is equivalent to roughly 21 pc in 
linear scale. NGC 4244 is nearly bulgeless, except for 
a compact and bright nuclear star cluster (see Seth et 
al. 2008 for a detailed study on the nuclear cluster; the 
bulge-to-disk ratio is 0.03 according to Fry et al. 1999). 
NGC 4244 exhibits systematic displacements of the mid- 
plane of a few tens of parsecs with a wave-length of a 
couple of kiloparsecs (this phenomenon has been termed 
'corrugation' and in the case of NGC 4244 it has been 
studied by Florido et al. 1991). Kodaira & Yamashita 
(1996) found NGC 4244 to have less warm dust and star 
formation than galaxies of the same Hubble type and H i 
content, thus labeling it as an 'anemic' galaxy. Finally, 
de Jong et al. (2007) found a sharp disk cut-off which 
they suggest to have a dynamical origin. NGC 4244 could 
be weakly interacting with NGC 4214, which has a simi- 
lar recession velocity and is found at a projected distance 
of 120 kpc from NGC 4244. 

Several studies have attempted to identify the struc- 
tural components of NGC 4244. Bergstrom et al. (1992) 
found no evidence for a massive stellar halo. Fry 
et al. (1999) detected only a thin disk. Hoopes et 
al. (1999) found that, unlike some other edge-on galax- 
ies, NGC 4244 has its ionized gas concentrated in the 
thin disk. Strickland et al. (2004) found no evidence for 
a radio, UV, or X-ray emitting halo. Seth et al. (2005b) 
found that the scale-heights of the older stellar compo- 
nents are larger than those of the younger components, 
which could be an indicator of the existence of a thick 
disk; however, the difference in scale-heights between 
the main sequence population and the red giant branch 
(RGB) population is smaller than in all the other galaxies 
studied in their paper. They also detected what could be 
interpreted as a halo. Tikhonov & Galazutdinova (2005) 
claimed to have discovered a tenuous extended halo with 
a scale-length of several kiloparsecs in NGC 4244. Seth 
et al. (2007) found a faint halo component 2.5 kpc above 
the mid plane at the position of the minor axis. 

To sum up, previous studies on NGC 4244 show little 
to no emission at any wavelength outside its well-known 
thin disk. As nearly all very inclined disk galaxies have 
been shown to host a thick disk, the absence of it in 
NGC 4244, if confirmed, would be an indication of a 
peculiar evolution or formation mechanism of this galaxy. 

In this Paper we present results of a study of NGC 4244 
based on a Spitzer Survey of Stellar Structure in Galaxies 
(S'^G) 3.6fim image (Sheth et al. 2010; Sect. 2), confirm- 
ing that this galaxy does in fact have signs of two disk 
components. For this purpose we have fitted the galaxy 
using Galfit (version 3.0; Peng et al. 2002; Peng et 
al. 2010; Sect. 3) and produced luminosity profiles of 
NGC 4244 at different galactocentric radii (Sect. 4). We 
describe the luminosity profiles and compare them with a 
simple inclined two-disk model with a diffuse dust contri- 
bution in the thin disk in Sect. 5. We present a summary 
of the NGC 4244 main structural components in Sect. 6, 
and discuss our findings in Sect. 7. 

2. PROCESSING OF THE S"'G 3.6AtM IMAGE OF NGC 4244 



The S'^G, described by Sheth et al. (2010), is a sur- 
vey using the 3.6 /xm and 4.5 /xm filters of the Spitzer 
Space Telescope Infrared Array Camera (IRAC, Fazio et 
al. 2004) with the aim of studying the stellar mass dis- 
tribution in a sample of 2331 local galaxies. 

The 3.6^m image of NGC 4244 we have used for our 
study has been processed by the S^'G pipelines, resulting 
in a mosaicked flat-fielded image and in a mask of fore- 
ground stars, background galaxies and bad pixels (Sheth 
et al. 2010). However, as we aim to reach the faintest pos- 
sible surface brightness level, we have manually refined 
the mask and closely examined the background. The 
background of the image is at a level of 0.0330 mJysr"^ 
(corresponding to ^3.6^m(^-B) = 24.2 mag arcsec"^). 
The cr of the background is of the order of 0.002 mJy sr~^, 
which is comparable to the offsets applied to the individ- 
ual frames in order to obtain a mosaic with a uniform 
background (0.0015 mJystr^^). We find a systematic 
background gradient which is shown in the top-left panel 
in Fig. [T] This gradient has an amplitude which is ap- 
proximately two times larger than that of the a of the 
background. The gradient is unlikely to be caused by zo- 
diacal light or by the concentration of bright stars at the 
northern part of the frame because it is not consistent 
with that of the S^G 4.5/im image which was obtained 
for the same galaxy. We have modeled the background 
using a third order two-dimensional polynomial fit to all 
unmasked pixels at a distance of 250 pixels (187.5") or 
more from the galaxy mid-plane. In this procedure there 
is no risk of contamination by any component of the 
galaxy, because we are going farther than any trace of 
extended component found by Tikhonov & Galazutdi- 
nova (2005) using the RGB star-counting method and 
because the surface brightness of the halo detected by 
Seth et al. (2007) drops below the S'^G limiting magni- 
tude at that distance from the mid-plane. We then sub- 
tracted the modeled background from the original image 
(the result is shown in the top- right panel in Fig. [T]). 

We rotated the image 43.6° counter-clockwise in order 
to place the mid-plane of the galaxy horizontally (lower 
panel in Fig. [ij. We estimate that the uncertainty in 
the rotation angle is about one degree. The rotation 
angle has been calculated by finding the mid-plane for 
all galactocentric distances r < 350arcsec (r < 7.5 kpc). 
This rotation angle is within the uncertainties consistent 
with the value of 42° used by Olhng (1996) and Fry et 
al. (1999). Once the galaxy has been rotated, we can 
define as the top part of the galaxy the area extending 
above its mid-plane and as the bottom part of the galaxy 
the region below it. This definition is arbitrary, but it 
will be useful because the galaxy is not symmetric with 
respect to its mid-plane. 

NGC 4244 has a warp-like feature starting from galac- 
tocentric radii of 4 to 5' in both optical (indicated with an 
arrow in Fig. [2 SDSS DR7 image; Abazajian et al. 2009) 
and near-infrared bands (indicated with arrows in the 
bottom panel in Fig.[T]) in agreement with the warp prop- 
erties measured by Saha et al. (2009). The color image 
of the galaxy (top panel in Fig. [2]) shows that this warp- 
like feature is delineated by a succession of blue knots 

knot ^knot) (^surroundings ^surroundings 

) ~ -0.4], 

which are areas of ongoing star formation as confirmed by 
the Ha image of Kcnnicutt et al. (2008). This indicates 
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Fig. 1.— Top panels: the 3.6/im S^G image of NGC 4244 before (left) and after (right) removing the gradient in the background. We 
have used the same histogram equalization display for both images. These images have been smoothed using a Gaussian function with 
a = 1.25 pixels. North is up and East is left. Lower panel: 3.6/.tm S^G image of NGC 4244 rotated 43.6° counter-clockwise in such a way 
that the galaxy plane is horizontal. In this frame North is roughly located at the top-left and East at the bottom-left. The image is 15' 
or 19.2 kpc in length. The distance between the nuclear cluster and the 'A' arrows is around 5'. The display for this panel is logarithmic. 
The arrow in the top-right panel indicates a light streak caused by the combination of the light of three foreground stars. The 'A' arrows 
in the bottom panel indicate the location of the warp-like feature. 



that the warp-hke feature is hkely to be a product of 
the presence of spiral arm fragments embedded in a very 
inchned galaxy disk (this possibility has been discussed 
by van der Kruit 1979). In the bottom panel of Fig. [2] 
the warp-like feature is delineated by the deviations of 
the tips of the gray areas (indicated with 'A' arrows) 
from the horizontal orientation. This warp-like feature 
is part of the corrugation pattern described by Florido 
et al. (1991), as outside the 'warp', at a lower surface 
brightness, the disk mid-plane bends in the other sense 
(seen in Fig. [2]) . In the top images of Fig. [T] (in which the 
color-scale distribution is based on a histogram equaliza- 
tion algorithm), the warp- like feature we have described 
may be confused with a light streak which goes downward 
from the south-eastern edge of NGC 4244 (opposite to 
the warp-like feature and indicated by an arrow). This 
streak is caused by three relatively bright aligned stars 



whose PSFs overlap in the histogram equalized display. 

The warp-like distortion we find coincides with that 
seen in the i?-band by Fry et al. (1999). These authors 
state that the optical 'warp' they see is coincident with 
the Hi warp discussed by Oiling (1996), but this conclu- 
sion is due to them having used an image of NGC 4244 
flipped along the mid-plane. The warp found in H i by 
Oiling (1996) starts at a galactocentric radius of r 10', 
which would place it at larger galactocentric radii than 
shown in the image in the bottom panel of Fig. [TJ 

An inspection of the SDSS DR7 image of NGC 4244 
(Fig. [21 indicated by the arrow labeled with letter 'B') 
shows a dust lane parallel to the mid-plane and slightly 
above it. This indicates that what we have defined to be 
the bottom part of the galaxy is facing us and thus we 
can describe it as the near side of the galaxy. This was 
already pointed out by Oiling (1996), who found that the 
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Fig. 2. — SDSS DR7 ugriz image of NGC 4244 oriented identically and at the same scale as in the bottom panel of Fig.[T] The top panel 
contains a color version of the image and the bottom panel has a black and white version of it. In the color image, the blue knots denote 
star formation and the reddish colors dusty areas. In the bottom panel, the 'A' arrows indicate the location of the warp-like feature which 
is described in Sect. 2. The 'B' arrow indicates a dust lane that runs parallel to the mid-plane of the galaxy and is located slightly above it. 



inclination of the galaxy is i = 84.5°. The dust lane is 
undetectable in the S^G images, even in a 3.6 fim—A.5 /im 
color-index map. 

3. GALFIT ANALYSIS OF NGC 4244 

In order to know whether a single-disk description is 
good enough to account for the structural components 
of NGC 4244 we have run Galfit (version 3.0; Peng et 
al. 2002; Peng et al. 2010) on the S-^G NGC 4244 im- 
age using the Galfidl procedures and display routines 
(see http://sun3.oulu.fi/~hsalo/galfidl.html). Galfit is 
a code designed for fitting one or more two-dimensional 
components to an image. 

As an input we used the background-subtracted S^G 
image, the refined mask, and a 'sigma' image, which ac- 
counts for the statistical uncertainties associated to each 
pixel. The 'sigma' image is made using Galfidl, which 
takes into account the pixel weights varying in different 
parts of the image mosaic. 

3.1. Inclined disk model 

In this model Galfit assumes a radially decaying ex- 
ponential profile. In order to run it we fixed the center 
of the galaxy at the position of the nuclear cluster and 
the disk PA as 46.4° which corresponds to 90° in the 
rotated images we present in Fig. [H We left as free pa- 
rameters the disk luminosity, the disk axis ratio, and the 
disk scale-length. As the galaxy appears to contain no 
bulge, no further components are a priori needed. The 
major caveat of modeling a very inclined galaxy this way 
is that the inclined disk Galfit model is an infinitely 
thin circular disk model, which will cause all the vertical 
extension of the galaxy to be interpreted as caused by 
the disk inclination. The resulting disk has an axis ratio 
d/D = 0.14, which gives a disk inclination of i = 82°. 
Considering that in our model the disk is infinitely thin, 
which will automatically cause the disk to be less in- 
clined than it really is, this value is compatible with the 
inclination oi i — 84.5° given by Oiling (1996). The in- 
clined disk model and the residual image (the original 
background-subtracted image from which the model has 
been subtracted) can be seen in the left column of Fig.[3l 



3.2. Edge-on disk model 

When applying the edge-on disk Galfit model we are 
assuming that the inclination angle of the galaxy is i — 
90°. The function which Galfit uses is 

where Eq is the central surface brightness, Ki is a Bessel 
function, r is the galactocentric distance, is the major 
axis disk scale-length, z is the distance above the plane, 
and hz is the scale-length perpendicular to the plane, as 
developed by van der Kruit & Searle (1981) for the light 
distribution in edge-on disks. As for the other model, 
we fixed the center of the galaxy at the position of the 
nuclear cluster and the disk PA as 90°. We left as free 
parameters the disk luminosity, the disk scale-length, and 
the disk scale-height. The major caveat of modeling a not 
fully edge-on galaxy this way is that we are ignoring the 
fact that part of the apparent scale-height is an effect of 
the galaxy inclination. The edge-on disk model and the 
residual image can be seen in the right column of Fig. |3l 

3.3. Comparison of the two NGC 4^44 Galfit models 

The main difference between the inclined and the edge- 
on model (see Fig. [3]) is that in the first case the galaxy 
model isophotes are elliptical and that in the second case 
they are diamond-shaped. The reason for this is that, in 
the first case, we are seeing the projection of an infinitely 
thin disk, which gives as a result a very eccentric ellipse. 
In the second case we are seeing a diamond shape because 
the Galfit light distribution for an edge-on disk has 
a cusp at the galactocentric distance r — 0. The real 
galaxy image is more peaked than the inclined disk model 
for galactocentric distances larger than r = 4' (which 
causes the butterfiy-shaped dark areas in the residual 
image) but has rounder isophotes than the edge-on model 
near r — 0. The residual image shows that, again, the 
inclined disk model is inaccurate at high galactocentric 
radii where it predicts a light distribution parallel to the 
z-axis much wider than observed. 

In both models, the residual image shows bright 
patches close to the mid-plane. In the residuals of the 
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Inclined single-disk model Edge-on single-disk model 
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Fig. 3. — The inclined single-disk model is presented on the left and the edge-on single-disk model is presented on the right. Prom top 
to bottom, the images correspond to the original S^G background-subtracted frame, to the Galfit model, and to the residual image. All 
the images are in a logarithmic display. The white rectangle indicates an arc-like feature that we interpret as possibly of tidal origin. 




Fig. 4. — Close-up of the box in the lower-right panel in Fig. [3] 
The proposed path of the arc-like feature is indicated with a black 
line. 



edge-on model, these patches are reminiscent of a small 
bar or/and a spiral arm system seen nearly edge-on. The 
residual image of the edge-on model also shows an en- 
velope of residual light which may be related to a thick 
disk. The extension of this envelope roughly corresponds 
to that of the galaxy's white-color envelope in the color 
panel in Fig. [21 However, part of this residual emis- 
sion, especially that at the top-right side of the image, 
has more to do with the warp-like feature described in 
Sect. 2 than with the presence of a thick disk, as in- 
dicated in Fig. [2] by the behavior of blue knots of star 
formation at high galactocentric radii. 

An arc-like feature is visible in the residual images of 
both models at the top-left side of the frame and at a very 
low surface brightness (Figs. [3] and HJ. It is more visi- 
ble in the edge-on model residual image (indicated with a 
white rectangle) and we speculate that it could be related 
to some kind of tidal feature, as its shape is somewhat 



similar to that observed in other very inclined galax- 
ies (such as in NGC 5907 studied by Shang et al. 1998 
and Martmez-Delgado et al. 2008 and in NGC 4013 by 
Martfnez-Delgado et al. 2009; these cases are admittedly 
much more spectacular than that described here). This 
feature is knotty, which can be explained by the presence 
of globular clusters or star-forming regions in a tidally 
disrupted galaxy. 

Neither of the models we ran reflect accurately the be- 
havior of the stellar components in NGC 4244. However, 
the simple modeling we have done is enough to show that 
a simple single-component description is not adequate to 
describe NGC 4244's structural properties. 

4. VERTICAL LUMINOSITY PROFILES 

In order to search for traces of a thick disk, we have 
produced luminosity profiles parallel to the z axis of the 
galaxy at different galactocentric radii (r). These lumi- 
nosity profiles are shown in Fig. [5l 

To produce the luminosity profiles, we first made ver- 
tical bins with a width of 225.0" (300 pixels or 4.8 kpc), 
and centered at galactocentric distances r = —337.5", 
-225.0", -112.5", 0.0", 112.5", 225.0", and 337.5". The 
225.0" bin width has been chosen because it offers a good 
compromise between spatial resolution and the profile's 
signal-to-noise. As the bin width is double the separation 
between bins, the bins are overlapping. For each bin we 
determined the surface brightness at a given z making 
a median of the counts excluding the pixels which were 
masked. We determined the location of the mid-plane of 
the galaxy from the profiles obtained in each bin. This 
was done by folding the profiles and minimizing the dif- 
ference between the upper and the lower halves of the 
light profile, with a precision of 0.1 pixel. For each bin 
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Fig. 5. — Luminosity profiles along the z axis at different galactocentric radii (seen in the subplot's titles). Short dashed lines indicate 
the fits to inner exponentials and the range of z at which the fit has been done (vertical lines). Dotted lines indicate the fits to outer 
exponentials and the range of z at which the fit has been done (vertical lines) for these profiles for which an outer exponential has been 
fitted (in some cases the vertical line indicating the inner boundary of the outer exponential fitting range is coincident with the outer 
boundary of the inner exponential fitting range). Dash-dotted lines indicate the minimum z at which the noise level has been calculated. 
The horizontal long-dashed line indicates the noise level. The profiles from the near side of the galaxy arc on the left and those from the 
far side are on the right hand side. The error bar length at the noise level is on the order of the line width. 



we made a visual guess of where the galaxy effectively 
ends. The noise level was calculated by producing the 
standard deviation of those pixels at higher z. This is a 
conservative estimate of the noise, as the profile is pro- 
duced by doing the median of the pixels in the bin. We 
smoothed the profiles in the vertical direction in order 
to get a higher signal at low surface-brightness areas. 
This smoothing has been done logarithmically and using 
a vertical scale-height of 5 pixels (3.75"). 

As the near and the far side of the galaxy are differ- 
ent, we have fitted them separately. We decided to make 
exponential fits at ranges in z where the data appear 
roughly exponential. The fit was made in magnitude 
units and all the points were given equal weight. Wc fit- 
ted one or two exponentials to each profile depending on 
whether we detected - by eye - any trace of a thick disk 
or not. The fitting ranges appear in Table [l] and they are 
indicated in Figs. [5] and |6l When preparing the fits, we 



considered all the luminosity values above the noise level 
to be reliable. The noise in our profiles is /i3.6^im(^S) ^ 
26.5 mag arcsec"^. If we had not taken the background 
gradient into account, as described in Sect. 2, this level 
would have been /X3.6;jm(^S) ^ 26.0 mag arcsec"^. The 
results of the fits are overlaid on the profiles in Fig. [5] 
and are summarized in Table [1] 

5. RESULTS 

We find that all but one of the far side profiles can be 
fitted with a single exponential. Only for the profile at 
—225.0" < r < 112.5 there is an extra faint component 
which cannot be fitted by an exponential. The near side 
profiles are well fitted by two exponential components 
except for that at r = 337.5" (Fig. [SJ. The scale- heights 
of all these fits are roughly stable with varying galacto- 
centric radius. 

5.1. Looking for a thick disk in the luminosity profiles 
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Figure \5\ (continued) . 
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The luminosity profiles (Fig. [S]) expose a lack of sym- 
metry between the far and the near sides of NGC 4244. 
Once an inner exponential has been fitted, it is clearly 
seen that larger z sections of the near side profiles can 
be fitted with a shallower exponential. By comparing 



lower and upper profiles, we can see that the concave 
infiection points between two exponentials in the near 
side profiles (which are especially noticeable for —112.5", 
112.5", 225.0", and 337.5") have no clear counterpart in 
the far side part of the galaxy. The infiection point is not 
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Fig. 6. — Representation of the fitting ranges in a 3.6/xm image of NGC 4244. The displayed area is the same as in the bottom panel in 
Fig. [l] The vertical solid lines denote the center of the bins. As in Fig. \5\ short dashed lines indicate the range of z at which the fit to the 
thin disk has been done. Dotted lines indicate the range of z at which the fit to the thick disk has been done. Dash-dotted lines indicate 
the minimum z at which the noise level has been calculated. The labels at the center of each bin indicate the bin galactocentric distance, 
r. 



SO evident in the near side profile at r = —337.5" and r ~ —225.0" and it is not noticeable at all for r = 337.5' 
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Notes: The uncertainties in this table are statistical errors (3a-). Galactocentric distance (col. 1), z above wrhich noise has been calculated 
(col. 2), noise level (col. 3), z range at which the fit for the inner exponential has been calculated (col. 4), inner exponential central surface 
brightness (col. 5), inner exponential scale-length (col. 6), z range at which the fit for the outer exponential has been calculated (col. 7), 
outer exponential central surface brightness (col. 8), and outer exponential scale-length (col. 9). 

o of a thick disk in the far side profiles, even at galacto- 

centric radii not affected by the non-exponential excess 
1^ .■ of light, indicates that a thick disk there is less obvious 

than in the near side of the galaxy, if it is present at all. 

The luminosity profiles show some anti-symmetry, 
which is expected for a more or less symmetric warp or 
a symmetric spiral arm pattern embedded in the thin 
disk. This anti-symmetry is seen in Fig. [SJ whose panels 
at large negative radius at the near side resemble most 
the panels at large positive radius at the far side, and 
vice- versa. 

The scale-height of both disk components is stable from 
r = -337.5" to r = 112.5" (Fig. E]). For this range of 
galactocentric radii, the thin disk scale-height is on av- 
erage 16.2" (350 pc) for the far side of the galaxy and 
14.1" (300 pc) for the near side of it. At the same range 
of galactocentric radii the thick disk has a scale-height of 
17.4" (370 pc). For comparison Seth et al. (2007) fitted 
a single disk with ^ 300 pc in scale-height. At galacto- 
centric radii larger than r — 112.5", both the thin and 
thick disks apparently flare (this can also be seen at the 
right side of Fig. [5] and has been described by van der 
Kruit & Searle 1981). This apparent flare is related to 
the warp-like feature, which we argued may be related to 
the spiral structure of the disk. 

The definition of scale- height used here is different from 
that used in other works such as the one by Yoachim & 
Dalcanton (2006) and thus the results are not directly 
comparable. In other papers a two-component model is 
fitted, but in our case, due to the very subtle change in 
the slope of the luminosity profile between the thin and 
thick disk-dominated z ranges, we have considered safer 
to fit a single exponential in each z-range. As a conse- 
quence, what we have defined to be the thin disk scale- 
height is less steep than what would be obtained from a 
two-component fit as it is influenced by the underlying 
thick disk component. 



Fig. 7. — Scale-height of the disks at different galactocentric radii. 
Dashed lines indicate the thin disk and the dotted line indicates the 
thick disk. Black lines indicate the measurements of scale-heights 
in the near side of the galaxy and red lines the measurements in 
the far side of it. 



may be due to some effect related to this region is being 
affected by the warp-like feature. The inflection point is 
also subtle at r = 0.0". 

An excess of light appears in the far side of the galaxy, 
for -225.0" <r< 112.5", but it cannot be easily fitted 
by an exponential and, even where we could in princi- 
ple have done it, the fit would have to be done over a 
much smaller range in z. This excess corresponds to the 
arc-like feature in the residual image obtained from Gal- 
fit (indicated with a white rectangle in the bottom-right 
panel in Fig. [3] and shown in Fig. IJ]). As it has a surface 
brightness comparable to that of the second exponential 
we find in the near side profiles, it could hide part of 
a thick disk in the far side of the galaxy for radii from 
r = -225.0" to r = 112.5". However, the lack of traces 



5.2. Thin and thick disks and opacities 
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Fig. 8. — Representation of luminosity profiles of a highly inclined 
(i = 84.5°) modeled galaxy with a semi-transparent thin disk with 
an /izT = 2.5/izt, loT = 0.5/ot, and a dust distribution in equilib- 
rium with the thin disk component with an opacity set in a way to 
have an optical depth of r = 1 at z = at the arbitrary galactocen- 
tric radius we are plotting, r has been measured through the line 
which goes from observer eye to the infinity beyond NGC 4244 in 
the direction we are observing. The far side profile is shown in red 
and the near side profile is shown in black. Dashed lines indicate 
the exponential fits to the thin disk and dotted lines indicate the 
exponential fits to the thick disk. The dash-dotted line indicates 
the light contribution of the thin disk and the long-dashed lines 
indicate the light contribution of the thick disk. The tick marks in 
the vertical axis are drawn every 0.5 magnitudes. The tick marks 
in the horizontal axis are separated by a distance equivalent to 
the fitted thick disk scale-length. 2 = has been set to be the 
mid-plane of the galaxy. 



We now explore whether the fact that NGC 4244 is 
not totally edge-on {i = 84.5° according to Oiling 1996), 
could explain the lack of traces of a thick disk in the 
far side of the galaxy. In other words, we want to test 
whether, except for the arc-like feature, the near and far 
sides of NGC 4244 could be intrinsically symmetric with 
respect to the mid-plane of the galaxy. To do so we have 
designed a model galaxy with two disks, each with the 
following light emission distribution per unit of volume 
(this light emission distribution can be considered as a 
mass distribution if we assume a constant mass-to-light 
ratio) : 



e(i?,z) 



(2) 



where eg is the central emission per unit of volume, R 
the distance to the rotation axis of the galaxy, z the 
height above its mid-plane, /ir the radial light emission 
distribution scale-length and the vertical light emis- 
sion distribution scale-height. The parameters of each 
disk will be from now on designated with subindices t 
for the thin disk and T for the thick disk, eo should not 
be confused with Iq, which is the fitted central surface 
brightness of the galaxy projected in the sky plane, and 
hz should not be confused with Hz-, which is the fitted 
scale- height of the disk in the projected galaxy. Even in 
the case of a galaxy with no dust and with i — 90°, Hz 
would be different to hz, as hz is similar to the definition 
of scale-height used by Yoachim & Dalcanton (2006) and 
Hz is the scale-height measured as a single-component 
exponential fit to a given range in z (see the discussion 



at the end of Sect. 5.1). 

The galaxy model code was designed to have a tunable 
thick disk central light emission, eoT, and a tunable thick 
disk mass distribution scale- height (/izt)- The disk scale- 
length was designed to be the same for the thin and the 
thick disk, which has been found to be true from obser- 
vations as a first approximation (the ratio between thick 
and thin disk scale- lengths is on average around 1.25 
according to Yoachim & Dalcanton 2006). We made the 
scale-height of the thin disk to be one tenth of the scale- 
length. We made the thin disk semi-transparent, with 
a local density of dust proportional to the local mass 
density, which means that the dust is assumed to be in 
equilibrium with the stellar thin disk component. This 
assumption is the most important caveat in our model- 
ing, as we do not know which is the real distribution of 
diffuse dust in NGC 4244. This imphes that the results 
in this section are only qualitative. We assume that the 
thick disk is transparent, i.e., contains little or no dust. 
The density of dust in the thin disk is a tunable parame- 
ter in our model. Our code projects the light emitted by 
our three-dimensional model galaxy using a tunable in- 
clination angle, i. We have studied our model at r equal 
to one scale-lenght of the disks. 

In order to get vertical luminosity profiles similar to 
that in Fig. [5] at the galactocentric radii at which a thick 
disk is detected, we have manually tuned the variable 
parameters of our model, except the galaxy inclination, 
which we have set to be z = 84.5° (Oiling 1996). We find 
that a good agreement is obtained when eoT = 0.5eot, 
hzT = 2.5/iztj and the dust density is set to have an 
optical depth r = 1 when z = at the arbitrary galac- 
tocentric radius in which we are making the luminosity 
profile (see Fig. [8|). The selected r is on the order of 
that deduced for the mid-plane of an edge-on galaxy in 
the iiT-band by Bianchi (2008) and thus is a reasonable 
upper limit for the mid-plane obscuration at 3.6/im. All 
T values in this paper have been measured through the 
line which goes from the observer to infinity, in the direc- 
tion of NGC 4244. Using these parameters, we find that 
for the bottom parts of the galaxy the ratio between the 
thick and the thin disk scale-height of the integrated dis- 
tribution of light along the line of sight is Hzt ~ 1.257?zt, 
and that the fitted central surface brightness of the thick 
disk is between 0.5 and 1.0 magnitudes dimmer than that 
of the thin disk, which is compatible with values obtained 
from the fits in Fig. \5\ 

In Fig. [Sj as in our model, the luminosity profiles gen- 
erally show a steeper thin disk scale-height for the near 
side of the galaxy (300 ± 15 pc against 350 ± 30 pc when 
the scale-height is averaged between r — —337.5" to 
r = 112.5"; the uncertainty corresponds to one standard 
deviation of the values used for the average scale-height 
computation). This is because for the far side of the 
galaxy, the thin disk opacity is obscuring the thick disk, 
thus removing light from the regions at low z, and mak- 
ing the slope in these regions (assigned to the thin disk 
in our fit) similar to that at higher z. This effect is seen 
in lines representing the thick disk luminosity profile in 
Fig.l 

In addition, for galactocentric radii between r = 
—225.0" and r = 112.5", the arc-like feature may be 
hiding any trace of a softened transition from the thin 
to the thick disk in the far side of the galaxy. At galac- 
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Fig. 9. — Representation of luminosity profiles of a highly inclined 
(i = 84.5°) modeled galaxy with a semi-transparent thin disk with 
/ij,T = 2.5/izt, loT = 0.5/ot, and a dust distribution in equilibrium 
with the thin disk component with an opacity set in a way to have 
an optical depth of t = 0.00, t = 0.50, t = 1.00, and t = 2.00 
(lines from up to down) at z = and at the arbitrary galactocentric 
radius we are plotting, r has been measured through the line which 
goes from observer eye to the infinity beyond NGC 4244 in the 
direction we are observing. The far side profile is shown in red 
and the near side profile is shown in black. The tick marks in the 
vertical axis are drawn every 0.5 magnitudes. The tick marks in 
the horizontal axis are separated by a distance equivalent to the 
thick disk scale-length obtained from the fit in Fig. [8] z = has 
been set to be the mid-plane of the galaxy. 



tocentric radii larger than r ~ 225.0", the galaxy disk 
starts to be distorted by, possibly, the presence of spiral 
arms, thus making the comparison between NGC 4244 
and the model more difficult. 

It can be argued that t = 1 is a very high opacity for 
a galaxy imaged in 3.6/xm, but this opacity only corre- 
sponds to the mid-plane of the galaxy, where the central 
dust lane is found, r drops fast with z, and it is on the 
order of r = 0.1 at the z at which the transition between 
the thin and the thick disk is found. In addition, we 
find that there is a degeneracy between the opacity and 
eor/eot- As an example, setting the mid-plane opacity 
as T = 0.6, and cqt = 0.65eot yields a luminosity profile 
very similar to that in Fig. El 

We have tested, in the case of a galaxy with i — 84.5°, 
Cot = 0.5eot and h^T = 2.5/izt the behavior of the lumi- 
nosity profiles as a function of the amount of dust in the 
disk. The result is shown in Fig. ^ in which we see that 
for a mid-plane optical depth of r = 0, both sides are 
symmetric. When increasing t, the central luminosity 
is reduced for both near and far side profiles. The re- 
duction in the central surface brightness causes Hzt (the 
measured scale-height of the thin disk) to be increased, 
making the change in slope between the thin and thin 
disk softer with increasing t. This effect is more pro- 
nounced for the far side profile. The trend described in 
this paragraph is general and also applies when studying 
highly-inclined disks with different eoT/eot and /izx/Zizt- 

We have also tested the effect of varying the galaxy in- 
clination, i. The results are displayed in Fig. [101 We see 
that for our modeled galaxy, in the cases of far side lumi- 
nosity profiles the effect of a thick disk is not enhanced 
due to inclination effects (the ratio of scale-heights be- 



z 

Fig. 10. — Representation of luminosity profiles of modeled galax- 
ies with difl'erent incUnations (i = 90.0°, i = 84.5°, and i = 79.0° 
which are represented in red, black, and green respectively) with 
a semi-transparent thin disk with an ft^x = 2.5hzt, ^OT = 0.5/ot, 
and a dust distribution in equilibrium with the thin disk component 
with an opacity set in a way to have an optical depth of r = 1.00 
when i = 84.5° at z = and at the arbitrary galactocentric radius 
we are plotting, t has been measured through the line which goes 
from observer eye to the infinity beyond NGC 4244 in the direction 
we are observing. The near side profiles are shown at left and the 
far side profiles are shown at right. Solid lines indicate the near 
side and dashed lines indicate the far side. The tick marks in the 
vertical axis are drawn every 0.5 magnitudes. The tick marks in 
the horizontal axis are separated by a distance equivalent to the 
thick disk scale-length obtained from the fit in Fig. [S] z = has 
been set to be the mid-plane of the galaxy. 

tween the thin and thick disk is constant with inclina- 
tion; for the model with aT = latz = and i — 84.5°, 
H^^/H^t = 1.20 for i = 90.0° and H^-r/H^t = 1-20 for 
i — 84.5°), but that in the case of near side luminosity 
profiles small angle differences with respect to the edge- 
on position make the difference between thin and thick 
disk scale-heights larger (for the model with a r = at 
z = and i = 84.5°, TJ^x/^^zt = 1-20 for i = 90.0° and 
HzT/Hzt = 1.27 for i ~ 84.5°). Larger angle differences 
with respect to the edge-on position makes any trace of 
thick disk to disappear. 

Our results in this Section confirm that the luminos- 
ity profiles of NGC 4244 are compatible with the galaxy 
having a thin and a thick disk and which are symmet- 
ric with respect to its mid-plane except for an arc-like 
feature on its far side. 

6. SUMMARY OF NGC 4244 FEATURES 
The features we have identified in NGC 4244 are: 

• An unresolved nuclear cluster studied by Seth et 
al. (2008). 

• Two disk components, whose identification can 
only be done in the near side of the galaxy, pos- 
sibly due to the presence of diffuse dust in the thin 
disk. 

• Possibly a spiral arm structure which causes the 
disk to have a slightly 'warped' appearance at 
around 5' from the nucleus. 

• Signs of a possible tidal feature. 
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NGC 4244 also has an H i disk which is coplanar with 
the optical and near-IR one except at very large galacto- 
centric radii (r > 10'; seen in Oiling 1996) where it shows 
a warp going in the opposite direction to the warp-like 
feature seen at ^ 4' — 5' in galactocentric distance. 

7. DISCUSSION AND CONCLUSIONS 

We find evidence of two disk components in a galaxy 
which was reported to have only one or as being a doubt- 
ful case. Using modeling we find that the fact of the thick 
disk is not detected in the far side of the galaxy can be 
explained by the effect of dust in a not perfectly edge-on 
galaxy 

Fry et al. (1999) did not detect the thick disk in their 
surface photometry because although the i?-band image 
has a limiting surface brightness similar to ours, they as- 
sumed that the galaxy was symmetric with respect to the 
mid-plane and they summed the top and bottom lumi- 
nosity profiles, which diluted the effect of the subtle thick 
disk. In addition i?-band is more sensitive to dust, mak- 
ing such a subtle thick disk hard to detect in this band. 
Seth et al. (2005b) deduced the presence of an extended 
component other than the thin disk from the scale-height 
of the oldest stars (RGB), which have a slight overden- 
sity above the fitted sech^^ profile. However, this com- 
ponent is more extended than the thick disk we found 
and may be a halo. The sum of all stellar components 
they have studied (young main sequence stars. Helium- 
burning stars, AGB stars, and RGB stars) shows some 
trace of what could be a thick disk in the online color ver- 
sion of their Fig. 3. Tikhonov & Galazutdinova (2005) 
found an extended feature by counting RGB stars at the 
location of the arc-like feature. Seth et al. (2007) stud- 
ied the RGB star distribution in a thin strip along the 
bottom part of the minor axis of NGC 4244 deep enough 
to find a halo component undetected by us due to its 
extreme faintness, but they did not notice any trace of 
the thick disk. This is probably due to the fact that the 
evidence for a thick disk is more subtle at r = 0" than at 
other galactocentric radii (see Fig. [5]). The scale-height 
for the disk they fit (~ 300 pc or ~ 14") is exactly the 
same as the one we find for the thin disk at r = 0". 

One possible thick disk formation mechanism is dy- 
namical heating of a thin disk. According to recent mod- 
eling (Bournaud et al. 2007; 2009), the thick disk should, 
at least in part, form due to gravitational instabilities. 
These are also responsible for the creation of giant star- 
forming clumps which will eventually merge to form the 
galaxy bulge during the first gigayear after the galaxy for- 
mation. As NGC 4244 has virtually no bulge and a thick 
disk which is not very differentiated from the thin disk, 
it is conceivable that for some reason the disk dynamical 
instabilities have been suppressed, which would lead to a 
lack of thin/thick disk differentiation. The other kind of 
dynamical heating is a consequence of the encounters of 



stars with disk substructure and is considered to be a sec- 
ular evolutionary process due to thin disk irregularities 
such as spiral arms and giant molecular clouds, or due 
to dark matter subhaloes. The lack of a very prominent 
thick disk would indicate that NGC 4244 has got a very 
smooth structure during most of its life. This smooth 
structure would have led to a modest 'blurring' in the 
disk of the galaxy. Alternatively, NGC 4244 could have 
been formed much more recently than other disk galaxies 
studied so far. 

Part of the thick disk could also have formed through 
accretion of stellar material from merging satellites at 
the time of the galaxy build-up (Yoachim & Dalcanton 
2006). According to this model the accreted gas builds 
the thin disk. In this formation scenario, the lack of an 
extended thick disk would imply that the merging proto- 
galaxies did not have a substantial stellar component at 
the time of the formation of NGC 4244. 

Last, but not least, the thick disk may have formed in 
situ. According to Elmegreen & Elmegreen (2006) edge- 
on galaxies with redshifts z = 0.5 — 4.5 (corresponding to 
look-back times of 4.8 — 12.2 Gyr) have disk morphologies 
whose ratio between long and short axes are compati- 
ble with those of local thick disks. However, these disks 
should shrink when the galaxy accretes gaseous material, 
which makes the potential well deeper. In the context of 
this thick disk origin theory NGC 4244 may be a case of 
a galaxy having accreted more gaseous material than av- 
erage, thus making the thick disk thinner than for most 
disk galaxies. This material accretion must have been 
through gentle gas infiows which have settled in the thin 
disk. 

To sum up, we have found evidence for two disk com- 
ponents in NGC 4244. The reason for the low degree of 
differentiation between the thin and thick disk remains 
unknown and may be related to a lack of strong secular 
evolution, to a lack of stellar material in the fragments 
which built up the galaxy, or to a high infiow rate of 
gaseous material after the galaxy formation. 
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